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1  Introduction 

1.  Historical  Background 

The  first  systematic  statistical  approach  to  the  problem  of  detection  in 
the  presence  of  noise  seems  to  be  that  of  Wiener  (Ref.  10),  who  studied  the 
problem  of  designing  a  linear  filter  for  separating  two  stationary  random 
processes  optimally,  using  the  least-squares  error  criterion.  Wiener's  work 
was  extended  by  Zadeh  (Kef.  11)  to  a  case  where  the  signal  wnich  it  is  desired 
to  filter  is  non-stationary.  Van  Vleck  and  Middleton  (Ref.  6)  also  treated 
the  problem  of  designing  a  "best"  linear  filter,  but  used  a  criterion  different 
from  Wiener's.  Their  work  was  getieralized  by  Dwork  (ref.  3).  These  analyses 
made  use  of  the  assumption  that  an  infinite  sample  of  observed  signal  is  availa 
ble,  and  adopt  a  more  or  less  artificial  criterion  for  the  optimum  detection 
process.  There  has  recently  been  scxne  work  on  analysis  of  finite  samples  of 
observed  signal  ^Ref.  2),  but  the  problem  of  obtaining  a  more  fundamentally 
acceptable  detection  criterion  has  received  only  scant  attention  in  the  engi¬ 
neering  literature. 

2.  Definition  of  the  Optirrum  Detector 

Before  defining  optimum  detector  it  seems  advisable  to  define  detector. 

Let 

S(t)  -  A  sin  (wt^e)  (1.2.1) 

be  the  signal  which  it  is  desired  to  detect.  A  and  u)  are  known  in  advance, 
and  9  is  purely  random.  (By  tnis  we  mean  that  9  varies  randomly  from  s/unple 
to  sample,  with  a  uniform  probability  density  in  the  interval  Let 

N^t)  denote  the  noise  voltage.  If  z(t)  denotes  the  observed  sample  of  length  T 
we  have 


z(t)  •  S(t)  ♦  N(t) 


(0  b  tf  T) 


(1.2.2) 
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Edgier  Reich  and  Peter  Swerling 


The  RAND  Corporation 
Santa  Monlca»  California 


*  ^  y 

This  paper  deals  with  the  problea  of  finding  the  ^optimua**  method  of 
detecting  a  sine  wave  of  known  frequency  and  amplitude  in  the  presence  of 
noise.  The  type  of  noise  considered  is  the  so-called  stationary  Gaussian 
process,  which  is  obtained  when  thermal  noise  is  passed  through  an  arbitrary 
linear  passive  device. 

The  analysis  takes  into  account  the  fact  that  in  practice  only  a  finite 
sample  of  observed  signal  is  available. 

The  optimum  detection  method  is  defined  as  that  which  maximizes  the  proba¬ 
bility  of  recognizing  the  presence  of  a  sine  wave  if  one  has  actually  appeared; 
while  the  probability  of  falsely  announcing  the  presence  of  a  sine  wave,  if 
none  has  actually  appeared,  does  not  exceed  some  pre-chosen  value. 

It  is  shown  that  when  the  noise  has  a  flat  spectrum,  all  the  relevant 
information  is  contained  in  the  amplitude  and  phase  of  the  Fourier  transform 
of  the  received  sample  at  the  frequency  of  the  sine  wave.  Almost  the  same 
result  holds  in  the  case  where  the  noise  has  an  exponentially  decaying  auto¬ 
correlation,  except  that  in  this  case  the  values  of  the  observed  sample  at  the 
endpoints  of  the  sample  also  play  a  role. 
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when  S(t)  is  present;  and 

t(t)  •  ^^(t)  (0<t  sT)  (1.2.3) 

when  S(t)  is  absent. 

Qjr  a  detector  is  meant  a  black  box  whose  input  is  z(t),  and  whose  output 
is  a  yes-no  indication  (Fig.  1),  indicating  the  detector's  guess  as  to  whether 
or  not  S(t)  was  present. 


leb 


~  no 


Fig.  1  Block  Diagram  of  Basic  Detector 

If  an  infinite  sample  were  available  (T  ■  oo )  it  would  be  possible,  in 
principle,  to  design  an  arbitrarily  good  detector  by  using  a  sufficiently  narrow 
band-pass  filter  centered  at  w ,  or  by  performing  some  such  process  as  looking 
for  a  periodic  component  in  the  autocorrelation  of  2(t).  Since,  however,  T  is 
finite,  it  will,  except  for  degenerate  cases  (sucn  as  a  missing  band  of  noise 
frequencies  around  w)  be  impossible  to  design  a  perfect  detector.  The  errors 
committed  by  the  detector  can  be  of  the  following  two  exhaustive  and  mutually 
exclusive  types: 

(a)  Detector  says  "yes"  even  though  S(t)  is  absent; 

(b)  Detector  says  "no"  even  though  S(t)  is  presentc 

> 

Let 

Pp  •  probability  that  the  detector  falsely  announces  the  presence  of 
S(t)  when  S(t)  is  not  actually  present. 

Pjj  •  probability  that  the  detector  recognizes  i(t)  when  3(t)  actually  is 


present. 


1 
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We  hav« 

Pp  -  probability  of  error  of  type  (a)  (1.2.4) 

Pjj  •  l-probability  of  error  of  type  (b)  (1.2.5) 

We  define  the  optimum  detector  as  the  one  which  maximizes  P|^  subject  to 
Pp  ^  pre-chosen  value  Pp*.  Generally  speaking,  the  higher  the  upper  limit  on 
is  set,  the  ligher  it  will  be  possible  to  make  Pp. 


3.  The  Likelihood  Ratio 

The  information  upon  which  the  detector  must  base  its  decision  is  all 
contained  in  the  observed  sample  z(t),  (O^tcT). 

Put 

-  k  J  (k  -  1,  2 . n)  (1.3.1) 

We  can  assume  first  that  the  detector  knows  only  ®nd  later  let 


n— 400.  It  is  physically  obvious  that  the  information  obtained  by  the  detector 
as  n — ^00  is  all  the  information  in  z(t).  This  has  been  shown  rigorously  in 
Ref.  4. 

In  the  general  case,  (z^,  z^,  ...,  z^)  is  a  set  of  random  variables  with 

a  certain  Joint  probability  density  which  depends  on  whether  or  not  S(t)  is 
present.  Let  U2,  ...,  u^)  -  probability  density  of  (z^^,  z^,  ...,  «^) 

when  S(t)  is  absent,  and  ^2’  '*n^  "  density  of  (z^,  Z2» 

...,  z^)  when  S(t)  is  present;  i.e., 

I 

L^(u^,  U2 . u^)  dUj^  du2  •  Joint  probability  that 

«1<  (i-l,2,...,n) 


when  S(t)  is  absent. 


(1.3.2) 
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L,  (u. ,  u^,  ....  u  )  du.  du.  . . .  du  >  Joint  probability  that 
1  1  <c  n  1  <  n 

“i  -  *i- “i^^'^i  .  *n) 

whan  5(t)  is  present,  (1.3.3) 

(L^  and  will  be  given  explicitly  in  section  II,  1.) 

The  detection  problem  is  equivalent  to  testing  the  hypothesis  that 
...,  z^)  has  density  function  against  the  hypothesis  that  it  has 
density  function  L^.  It  ft  a  known  theorem  in  the  theory  of  testing  statis¬ 
tical  hypotheses  (Ref's.  1,  U,  5»  6,  9)  that  in  order  to  maximize  subject 

to  Ppi  Pp*  one  proceeds  as  follows: 

(a)  Define  the  "likelihood  ratio" 


h^*l'  *2'  '*•'  *n^ 

A(*i.  *2 . *n^  "  - : 

^0^*1’  *2'  ’**'  *n 


(1.3.4) 


A  is  a  function  of  the  random  variables  z^^,  Z2$  ...»  and  thus  is 


itself  a  random  variable. 


(b)  Let 

•  probability  density  of  A  under  the 

hypothesis  that  S(t)  is  absent.  (1.3.5) 
W^(x)  •  probability  density  of  A  under  the 

hypothesis  that  S(t)  is  present.  (1.3.6) 
Knowing  the  functions  L^(Uj^,  U2,  ...»  u^)  and  Lj^(Uj^,  U2 . u^),  and 

using  (1.3.^.  and  can  be  found. 


(c)  Find  u  •  u(Pp*)  such  that 

00  ^ 

/  W^(>.)  dX  -  1  -  /  w  (X)dX  -  Pp*  (1.3.7) 

u  0 
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(d)  For  «  given  observed  sanple  *^)»  l«t  the  detector  s«j 


•yes”  when 


A(z^i  Z2*  ^ 


(I.3.B) 


and  "no"  when 


^2*  •*•»  ^ ^ 


(1.3.9) 


Thus,  by  (1.3.7),  the  probability  of  falsely  announcing  a  signal  to 
be  present  will  be  exactly  Pj.*, 


The  resultant  value  of  P|^  ^11 


Pj^-/Q,(X)  dX 


(1.3.10) 


V 


II  PTivation  of  Optimum  Detector 

1.  Noise  with  Arbitrary  Autocorrelation 

The  noise  N(t)  is  assumed  to  be  a  stationary  Gaussian  random  process 

2 

with  mean  0;  variance  a  ;  and  autocorrelation  function 


<  (r)  -  E[N(t)  NCt4t)j  -  E|N(t)  N(t-r)j 
If  signal  is  present,  it  has  the  form 


(1) 


(II. 1.1) 


S(t)  •  A  sin  (0t*9)  (II. 1.2) 

where  A  and  u>  are  known  constants  while  9  is  equally  likely  to  be  anywhere  in 
the  interval  0  <  0  s  2n. 

The  total  voltage  is 

z(t)  •  N(t) 

7.(t)  -  N(t)  ♦  S(t) 

iT 

Let  this  voltage  be  observed  at  times  “ 

and  denote  z(t^)  by  z^. 

We  will  assume  that  the  matrix "  l|2,...,n) 

is  non-singular  (i.e. ,  has  an  inverse)  for  each  n. 

If  signal  is  present,  the  conditional  probability  density  for  (Zj^,  ...»  z^) 

when  the  phase  is  9  is  given  by'  ' 


(signal  absent)  (II. 1.3) 

(signal  present)  (II.l.Z^) 

(i"l,2,...,n) 


•  •  •  • 


(2n)2l^ 


exp  < 


"  I  sin^u)tj^^e)j  [uj-A  sin(fc)tj«9)  f  (11,1,5) 


means : 


expected  value  of  X. 


Sometimes  this  is  denoted  by  X, 
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wh«r« 


[^j]  -[V] 

Kl-  determinant 

[!J  •[?.)•«■'■ 


Thus  the  Joint  density  function  for  Z2,  ...»  when  signal  is  present 


is 


2n 


u  d« 
n' 


(II.1.6) 


If  signal  is  absent  then  the  Joint  dsnsity  function  for  (z^,  Z21  ...»  z^) 


is 


L^(Ui,U2,.. 


n 


1 _ 

n  1 


(2Tt)2|f|2 


ij“i 


(II. 1.7) 


Hence  the  likelihood  ratio  A  is 

n 


I^(zj^yZ2»...»  Zjj ) 

^o^*l**2*  * ’ ' *  *n^ 


(II. 1.8) 


or 


A  • 

n 


1_ 

2n 


sinv(ot 


1*"’]  h- 


A  sin^u^t 


'ij 


or 


-b  2n 


■  - -  /  exp<  c  sin  0  ■»  d  cos  ^  ♦  g  cos  2©  -  h  sin  29  id©  (II.  1.9) 

n  2tt  y  [  n  n  ''n  n  J 
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where 


2  n 


’n'T  co.a;(t^-tj) 


'n  ■  *  *1 


n 


flinwtj 


(II. 1.10) 


,2  n 


*n‘- 


2  n 


’n-T 


The  integral  appearing  in  (II. 1.9)  is  cleerly  a  continuous  function  of 

hfCadaff  g_  . 

n’  n’  n'  n’  ®n 


Assuae  that,  as  n— >oo,  we  have 


b  — 
n 


c  — >c 
n 


(3) 


d  — >d 
n 


«n— 


h — ►h 
n 


(II. 1,11) 


(3) 

'  Finite  limits  b,  c,  d,  g,  and  h  do  not  exist  for  all  possible  autocorrelation 
functions.  In  most  interesting  cases,  however,  the  lirrits  will  exist.  (Since  c 
and  d  are  functions  of  the  random  variables  z^,  we  can  only  Sriy  that  c  and  d  ex?st 
with  probability  one.)  The  existence  of  these  (finite)  llniiis  is  closely  connected 
with  the  possibility  of  perfect  detection;  if  these  limits  do  exist,  perfect  detec¬ 
tion  is  impossible. 
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We  will  then  have 


» • • »  » )  — >  A 
n  1  n 


|ta)|  -  M(c,d;  b,g, 


h^ 


(II. 1.12) 


where 


-b 

M(c,d;  b,g,h)  -  ^ 


2n 


{ 


exp  <c  ain  ft  ♦  d  cos  6  ♦  g  coa  2ft  -  h  ai 


n  2ft| 


dft  (II. 1. 


Here  c  imd  d  are  fVinctionala  of  the  obaerved  z(t)  and  are  thua  random 
variablea. 

If  aignal  ia  actually  present,  the  vector  (c,d)  has  a  density  function  which 
we  shall  denote  by  F^^u,v);  if  signal  ia  actually  absent  (c,d)  has  a  density  func* 
tion  denoted  by  F^(u,v), 

According  to  what  was  said  in  Part  I,  the  optimum  method  of  detection  ia  as 
follows: 

(a)  Find  F^(u,v).  Thus  find  ^via  II. 1.12  and  11.1,13)  the  probability 
density  for  A  when  signal  is  absent. 

(b)  Let  Pp*  •  pre-chosen  upper  bound  on  Pp. 

Choose  ^(Pp*)  such  that  Prob  |^u^A  [signal  absentj  ■  Pp*. 

(c)  For  a  given  observed  sample  z(t\  compute  A  and 
U)  if  A  s  u(Pp*)  say*  signal  is  not  present. 

(ii)  if  A>  M(Pp*)  say  signal  ^  present. 

The  intervals  0<A<4  correspond  to  regions  in  the  (c,d)-plane  which 
can  be  determined  from  11. 1.13.  Thus,  the  above  procedure  is  entirely  equivalent 


13  > 


to  the  following  one: 
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Find  (numerically  or  otherwise)  the  regions  ri^  in  the  (c,d)-plane 
corresponding  to  0<  A  <  u. 

(b)  Find  F^(u,y)  -  the  dmnsit.v  function  for  (c,d)  when  signal  is  absent, 

(c)  Choose  R  so  that 


F^^u.v)  du  dv  -  Pp* 


jkcknotes  the  whole  ^c,d)-plane; 


L' 


R-R 


R-R  ■  complenjent  of  R  . 

J 


(d)  For  a  given  observed  sample  z(t),  conpute  c  and  d  and 
(i)  if  (c,dW  R^  say  signal  is  not  present. 


(ii)  if  (c,d^^  R^  say  signal  ^  present. 


The  first  three  of  these  steps  do  not  require  any  observations  to  be  made; 
they  can  be  carried  out  once  and  for  all,  once  the  autocorrelation  function  of 
noise  is  known. 

A  fifth  step  which  would  be  of  interest  is  the  calculation  of  the  probability 
of  detection  Pj^  (which  depends  on  Fp*-), 


(e)  Find  F^(u,v)  -  the  density  function  for  (c,d)  if  signal  is  present. 


Then 


r 

Pj^iPp*-)  -  j  F^(u,v)  du  dv. 


h-U 


The  calculation  of  F^^u,v)  and  Fj^(u,v); 

n 

E 

1 


We  have  c 


n  “  E^ij^i . j 


coaatt 


n 


1 


n  ‘-r  ^ij^i  sinwtj 
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(tt)  Signal  not  present  Calculation  of 


the 


2 

are  normal;  mean  0;  Tariance  o  ■  ^(0) 


c  ,  d  are  normal  with  mean  0*  Also 
n  n 


^  t.i.. 


2  2 

•  A  ^ - -  T.  4  .  2. 1,  cos  tut,  cosKt. 

"  l.J.K./  iJ  i  J  ^ 


^  #  J  •  *  »  • 


^  #  j  ^ 


E  rj^  co.Wj  co,-t^  (,lnc.[fjjJ.U  r^) 
j  •  ^ 


Thus 


-  A^  ^  T^j  coswt^  cosiutj  •  2{b^*g^)  (11,1.14) 


similarly 


— o  o  " 

d  -A  IZ  ,  sinu»t.  sinu)t,  •  2(b  -g  )  (11.1,15) 

n  1  ^  n  n 


and 


0  d  -  a 
n  n 


cos  wt,  sin  u)t4  •  2h 

X  j  rt 


(II. 1.16) 


Hence,  in  the  limit,  (c,d)  has  the  density  function 


•  ^7^  '*p 


1  r  2  -  2 

-  2  :“ll'^  *^^12  ^^^“22  ^ 


where 


and 


■  determinant 

2ib^g) 

2h 

2h 

2(b-g)_ 

®11  ®12 

?(b*g) 

2h 

®21  ®22 

2h 

» 

2(b-g) 

m 

-1 


(II. 1.17) 
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(b)  Signal  is  present  Calculation  of  Fj^(u,v^ 

2 

The  are  normal  with  means  A  8in(6<t^*ft);  variance  o  •  -tCO). 
For  a  fixed  9  we  now  have 


c  •  c  (9)  ■  ^  Xi  4  sln^Ort,  *9)  coe  UJt  (II. 1.18) 

n  n  ,  uj  1  j 


i  -  d  (9)  -  A^  Z2t4«  8in(U/t.»9)  slnwt, 
n  n  1  ^  ^ 


(II. 1.19) 


c  -c  "A 
n  n 


n  r  1 

C^ij|»i*'^  8in(U)ti^9^J  cos  Wtj  (II. 1.20) 


d-d  "a  53^4  *4-h  sin^UJt ,♦«)  sinWt.  (II. 1.21) 

n  n  ,  ij  1  1  •!  j 


But  the  variables  z^-A  sini<Vt^*9)  have  the  sare  statistical  properties 

as  have  the  wheifeonly  noise  is  present 

Hence  (c-c)  ,  (d-d)  ,  and  (c-c)(d-d)  have  the  same  values  as  when 
signal  is  absent.  The  only  difference  is  in  the  means: 


-  2 

:  {9)  •  k  8in(wt.-»9)  cosuJt, 

n  ^  aj  1  J 

2  p  ”1 

■  a  5r  ^  ^  cos  0)  tj  sin  Oit^  co89-»co8  (c)  t  J  sin9 


c  (9)  -  2h  cosO  ♦  2(b  ♦g  )  sin  9 
n  n  n  n 


Henc«, 


similarly, 


c(<i)  ■  2h  co8©*2(b^g)  8in0 


(1(0)  ■  2h  sin0  ♦  2(b-g)  co80 
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(IU.22) 


(11.1.23) 


The  conditional  density  liinction  for  (c,d)  when  the  phase  is  0  is  then 
Fj^'(u,v|0)  -  - ^  exp(-  [u-c(«)j  [u-d(»)]  ^1“  (11.1,24) 


where  D,  a^j  have  the  same  values  as  in  II.  1.17. 

Methods  for  explicitly  calculating  functionals  of  the  type  c  and  d  for 
gene  ral  autocorrelation^  srs  outlined  in  Appendix  II.  Sometimes  it  is  possible 
to  calculate  b,  c,  d,  g,  h  by  explicit  calculation  of  b^,  c^,  d^,  g^,  h^; 

an  example  of  this  appears  in  II. 2.  Once  c(0)  and  d(0)  have  been  obtained, 
th.e  density  function  F^(u,v)  for  (c,d)  in  the  presence  of  signal  is  given 


fi' 


(11.1.25) 


Fj(u.t) 


II.  2.  Noise  with  Exponentially  Decaying  Autocorrelation 


In  this  section  we  consider  the  case 


<  (c)  -  e 


-a  |c| 


(a,  ^^>0) 


(II. 2.1) 


^e  will  obtain  expressions  for  b,  c,  d,  g,  h. 

It  is  possiole  to  calculate  explicitly.  We  have 


ri 

L“ 


Y 

1 


-Off 


y 


^n-1 

p-2 


(II. 2.2) 


.  ^  -cd  »  T 

where  V-  e  ;  d  ■  - 

n 

As  can  easily  be  verified,  the  inverse  matrix  is 

1  0  0  ... 


(1-/) 


-/  1 


♦  0  ... 


n 

0 


•3  l»^  -y  ...  0 


0 

0 

n 


0  0  0 
0  0  0 


0  ...  -< 
0  ...  -y  1 


(II. 2.3) 


where  the  main  diagonal  is  1,  I*  ^ ,  1*3^,  ...»  1*3^»  1  ;  th^  diagonals  inurediately 
above  and  below  the  main  diagonal  consist  entirely  of  and  all  other  elements 

are  0. 

The  functionals  b  ,  c  ,  d  ,  g  ,  h  (II. 1.10),  the  limiting  values  of  which 

n  n  n’  °n’  n 

must  be  calculated,  all  contain  a  factor  of  the  fora 


n 


*  f  ^IJ  *1  yj 


,11.2.4) 


where  Xj^  ■  xCkJ),  yj^  -  ylki);  '^nd  where  at  least  one  of  the  functions,  soy 

y(t),  is  a  bounded  analytic  function  of  t  while  t.oe  other  function,  say  x(t), 
is  a  bounded  continuous  function  of  t. 


Putting  ^2o3)  into  (2.4)  t^ives 


B  - 

n 


I 

Using  the  relations 


2  ^ 

Vi^Vn’^^'^r)  E  Vk-  ‘  F  Vk.i-‘ t  V-i': 


-1 

r 

1 


d-y  A 

2  ^k  1 

c 

♦ 

.  je 

r'  ko 

A 

2  ^k  1 

where  v  - 


t  •  k^ 


,  etc., 


we  can  writ* 


n-1 

•1  •  ^  Vk 


?  Vk 


d  V  •■ 

2  Y  Vk 


..  ^2 

V  ♦  *— 

•'k  t 


n-1 


s'  ^  i’ V  *  ^  ^  tf  ^ 

^>k.i^  ■  Vk  -  Vk  *  T  "I'Vk  -  r  f  v«  • 


Putting  (2.6),  (2.7)  into  \2.5)  and  collecting  terms  gives 


S:  I^Vk  ?Vk  •^‘‘-’''^Vk'  Vn' 


pa->‘) 


r 

■’<^‘Vn-  *1^1'  •  2  '<n'n'  Vl’  '  c  'VV  Vl 


This  can  be  rewritten  as 


“n  •  fh  ?Vk 


2j7^^Vk  •  S  ‘Vr  Vn’ 


•:^  vx  -  X  ^- )  ♦  r(-^ 
2np  n  n  rl  n 
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(II. 2.5) 


0(^^)  (II. 2.6) 

(11. 2.7) 

(11. 2.8) 

.  r^y3) 


(II. 2.9) 
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bo  if  b  •  llm  b  , 
n  ' 

n— KD 


i  j 


y  (t)  dt  -  xtO)  y'(n)-ay(0)1 

“  J  °  L  J 


ai.2.10) 


♦  xtT^[y'(T)  ♦  aya) 


We  are  now  in  a  -csition  to  evaluate  the  liTiitintj;  valuea  of  the  func¬ 


tionals  ^II.l.ll): 


.2  r  >  2 

,  A  a/,  u;\_l 

b-r-  '^Cl*”::}!^- 
U  2p '  2 '  0 

0 


^(l  -»  /  *(t)  coswt  dt  ♦  ^  2^T)  fa  COSU. 


*•  / 

a  7. 


r  - w  9 i nwT 


^■576 


2 

'  z(t)  sinu^  dt  -  z\.C)  *  »(T)r4fc;«un,  ♦  0  sir. 


2/ 

a  J 


(iio2ai) 


il-OL/i 


2  1 

sin  2tfT  ♦  2iihr)j 

“  J 


V 
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II„  3.  White  Noise  ^with  Flat  Spectrum) 

A  definition  is  in  oroer  here  since,  strictly  speaking,  there  is  no  sta¬ 
tionary  Gaussian  random  process  with  a  spectrum  that  is  flat  for  frequencies 
from  0  to  ♦  ooo  Still  it  is  possible  to  "approach”  white  noise,  in  the  sense 
that  it  is  possible  to  define  a  sequence  “  1»2,...)  of  Gaussian  random 

processes  with  the  property  that  their  power  spectra  approach  a  flat  spectrum. 
Thus,  consider  white  noise  with  the  power  density  spectrum 

G(f)  ■  K  volt^/cps  { '  5  f <00  )  (11,3.1) 

This  will  be  defined  is  the  "limiting''  case,  as  n->oo,  of  a  sequence  (v‘! 

of  Gaussian  random  processes  with  autocorrelations  ~  e  .  ^Set 

Kef„  5,  pp  42,  for  relation  between  power  spec’ rum  and  autocorrelation.)  Putting 

a  ■  m,  d  “  V  (11.2,11)  and  letting  m-^oo  gives  the  following  limiting  values: 
U 

h 

^  "  Tk 


^  f  z^t)  cos<*^  dt 

^  J 

0 


^  t  z(t)  slnan,  dt 


(II. 3. 2) 


« • 


h  -  7—7  (1  -  cos  2iJl) 
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Note  thiit  the  endpoints  of  the  observed  sample  no  longer  play  a  part  in  the 
optimum  detection  px*oce88, 

A  particularly  simple  result  is  obtained  when  T  is  an  integral  multiple 
of  the  half  period,  of  S(t).  The  remainder  of  this  section  is  devoted  to 
this  case.  Note  that  g  and  h  /anish.  Thus  the  likelihood  ratio  becomes  (11.1.12^ 


A 


(II. 3. 3) 


where  I  is  the  modified  Bessel  function  of  the  first  kind,  of  order  zero, 
o 

2  2 

From  ^11.3.3  '  we  see  that  the  only  significant  statistic  is  ■*d  \  which 
is  proportional  to  the  s  .uared  amplitude  of  the  Fourier  transform  of  the  observed 
sample  at  radian  fre  iuency4J  An  explicit  formula  for  the  bias  level  ^  can  be 
obtained  for  this  case.  Since  1^  is  a  monotonic  function  of  its  argument,  the 

critical  value  for  A  corr<*spond3  to  a  critical  v-ilue,  s.y,-i3^Pp*)  for 

2  2 

V  •  (c  ♦d  ).  Wiien  5(t)  is  absent,  the  orobability  that  V  is  between  v  and 
v^dv  is  iRef.  1,  pp  236) 

v 

~  e  ,  tils  defined  as  a  function  of  P,>  bv  the 
40  r 

relation  correspondin.  to  3.7): 

(S  v 

dv  -  Fp*  (11,3.4) 


or 


(II. 3. 5) 
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Summing  up,  the  rule  for  this  optimum  detector  in  the  case  of  white  noise, 


n 


when  T  an  integral  multiple  of  ^  ,  is: 


z(t)  e'  dt 

o 


Announce  absence  of  S(t)  if 


2^t)  e  “  dt 


a_'\ 

.y) 


(II. 3. 6) 


2  2 

To  calculate  Pj^  it  is  necessary  to  obtain  the  probability  density  of  c  ♦d 

when  S(t!'  is  present.  The  conditional  probability  density  for  phase  0  turns 
out  to  be  independent  of  0;  it  is  equal  to 

V 


Prob  Vv£C^«d^ 


<v»dv|  •  e”^  e 

where  b  is  given  by  (H^3.2).  Therefore, 


(11.3.7) 


J  e”^  dv 


(II. 3. 8) 


-InPp* 


For  small  ratios  of  average  received  signal  energy  to  noise  poweiy^cpa, 
we  have  C<b<<l  ;  UI.3.8)  then  intnlies 


A^T 


(II. 3. 9) 


In  general,  it  follows  from  ^11.1,8)  that  if  a  * 


V^F 


(II. 3. 10) 


Of  course,  ^‘d  H-IO)  hold  for  .my  noise  autocorrelation  function  and 

any  value  of  T  if  the  detector  is  desirned  on  the  basis  of  the  theory  of  testing 
statistical  hypotheses. 


P-305 

f -11-52 

-21- 


hPPENDIX  I 

The  foregoirv’  can  be  formally  generalized  as  follows:  f^uppose  the  signal 
to  be  detected  is  ‘^Iven  by  3(t;aj^, . . .  -  3(t;a?)  where  are  K  parameters 

having  probability  districution  function  0(?'.  Let  the  noise  be  stationary, 
Gaussian  with  autocorrelation  1  (V)^  Then 


. ■  — n/2  --'irr 


I  iAi.n 


Therefore, 


i  i  ^  n  *- 1\(  ^ 


St  /u  -J(i. 


y-iG,^) 


Also  L  .  ,u  '' 

0  1  ’  n 


(2n^ 


n/2,,  |l/2  *i 


(al 


1  ^ 


1  n  ^ 

exp^-  2  j  J 


-  zr 


f i  /  \  V  Al .  i*  ' 


This  can  be  applied  to  the  case  where  are  all  fixed  instead  of  being 

randoBv U . ®a »  signal  shape  md  phase  are  known). 

Then  (denoting  the  signal  simply  by  G(t)) 


'n  ■  "lb 


■  «p  - 1 


n  n 

-AtJ  3(t,'/»  »‘Kp- -  Slj  Gvt^)  Zj 


(AI„5) 


Thus 


the  statistic  which  it  Is  nec'-ss-ry  to  calculate  i-*  St"j  Htp  Zj. 
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If  the  noise  is  white  noise,  this  (Quantity,  as  n— >00  ,  approaches  a  quantitj 
proportional  to 

T 

^  f  Cvt)  z(t)  dto  In  other  words,  all  the  information  ae  to 

^  J 
0 

whether  or  not  the  signal  is  present  is  contained  in  the  cross-correlation 
between  sirnal  and  observed  sample. 
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ArrsNDix  II 

The  problem  of  evaluatin^^  the  quantities  b,  c,  d,  r,  and  h  can  be  formally 
reduced  to  the  problem  of  solving  \  cert:  In  inte :ral  equation. 

Let  the  correlation  function  of  noise  be  given  by 


L‘\sa)  -f^9,t)-fit,5) 

The  quadratic  fonts  to  be  evaluated  are  of  t  form 


a' 


iAiI.l) 


n  •  Yl 


^AII.2) 


where 


f  it  t  J 

-  ^  J  - 


-I 


Xj  •  xitj',  where  x'.t''  is  a  bounded  continu  us  function 
in  0<t<T 

y^  •  y(t^  ,  where  y(t)  is  a  function  possessing  derivatives 
of  all  orders. 

Un  some  cases  we  may  have  yit)  «  x^t''.  ) 

Now  let 


then 


n 

S  ^  y  -  w 

“ij  j 


Zf  ,  -  y: 


^AII.3) 


r 


n 


Hence  we  must  evaluate  B  •  w.x,  where 

"  j-i  '  J 


{■] 


^AIIo3  / 


I  is  tnf  soluli  )n  of  1,AII.3K 


As  n — ^00  iand  as  tne  maxlmimi  interval  t^  t 't  ^  the  rroblem  formally 

reduces  to  the  evaluation  of 


TZT 


In  the  stationary  case  'fvS,t'' 


>s-t ). 


T*0  , 

/’  (5) 

/  x(t)  dW(t) 
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IAII.4) 


wri®re  W^t)  is  the  solution  of 


yit'  •  f  vs.t''  dVs’i,s)  V  <t<T^ 
•Ln 


(All. 5) 


As  an  excunple  consider  the  Ctise  <'(3,t)  ■  pe  discussed  in  section 

11.2.  L^-t  us  eviiluate  c,  for  exanple.  In  t'.is  case  y^t''  -  A  cosuA;  x(t)  ■ 

/.  ,t}.  (/illc5/  becomes 


T*0  ,  ^ 

/  .  /'  "0  h-t 

rt  COS^t."  J  rie  I 
-0 


dw  (s ) 


(All. 6) 


Tne  solution  W^t''  is 


2 .  t 


Wvt'  •  |t[i*  ^  J  COSkB  ds  ♦  R(t) 
''  L  a  0 


where  hit)  is  a 


function  which  is  ever'.’where  constant 


except  for  a  jump  of  it  t-O  and  a  Jump  of  [cosu/T  -  -  sinu/T  at  t"T, 

2p  a  _ 


Hence 


Ttn 

y  Zvt)  dW(t) 

-0 


A  A 

•:rr  — -  z(T)  a  cosa/r-<c;slnuT 

i.flr' 


2  T 

♦  /  zvt)  cos(ot  dt 

Q  J 

0 

which  agrees  with  the  value  previously  obtained  i3ee  II.  2. 11), 


means  lim  / 

V  t— 

— O  -( 


To  within  an  arbitrary  additive  constant 
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Our  rtblllty  to  ?olve  t be  inte»’r.<l  e-.u^ti-in  in  closed  fom  was 

8<}rr.ewh&t  fortuitous.  The  luthors  are  investigating;^  the  corresponding  equation 

-in 

for  autocorrelation  functions  of  type  ^  vt''  ■  fTp  e  ,  where  P  is  a  roly- 
nomial.  The  possibility  that  general  autocorrelation  functions  might  be 
approximated  by  Laguerre  polynomials,  and  npt roximations  to  the  values  of  c, 
d,  etc.,  obtained  in  this  way,  is  being  investigated. 


bh/P3 :  ime 
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LIST  OF  STKaOLS 


S(t) 

Nvt) 

z(t) 


V 

V' 

f  vD 
1 


ij 


A  sinvwt»«)  ■  sine  wave  signal 
noise 

observed  signal 

duration  of  observed  sample 

probability  that  detector  falsely  announces  the 
presence  of  S(,t'i  when  S^t)  is  not  actually  pr-sent, 

upper  limit  set  on  Pp 

probability  that  detector  recognires  S(t''  when  S(t) 
actually  is  present. 

^  ik  '  1.2,.  .,n) 

likelihood  ratio  with  the  Values  zit^^,  ii  *  l,...,n) 
as  the  observed  variableSo 

lim  A 

n 

n->x 

critical  value  of  A  •  bias  level 
a  r.onotonic  functi-^n  of  A 
critical  value  of  V  ■  bios  level 
autocorrelation  of  N(,t' 


h  ,c  ,d  ,g  ,h 
n  n  n’'^n  n 


functi  inals  occurring  in  calculation  of  A 


n 


b,c,d,g,h  •  limits  of  above  as  n->a> 


VjJ 
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